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Abstract 

A quantum field theoretical approach to the thermodynamics of dense Fermi 
systems is developed for the description of the formation and dissolution of quan- 
tum condensates and bound states in dependence of temperature and density. 
As a model system we study the chiral and superconducting phase transitions in 
two-flavor quark matter within the NJL model and their interrelation with the 
formation of quark-antiquark and diquark bound states. The phase diagram of 
quark matter is evaluated as a function of the diquark coupling strength and a 
coexistence region of chiral symmetry breaking and color superconductivity is 
obtained at very strong coupling. The crossover between Bose-Einstein conden- 
sation (EEC) of diquark bound states and condensation of diquark resonances 
(Cooper pairs) in the continuum (BCS) is discussed as a Mott effect. This effect 
consists in the transition of bound states into the continuum of scattering states 
under the influence of compression and heating. We explain the physics of the 
Mott transition with special emphasis on role of the Pauli principle for the case 
of the pion in quark matter. 

PACS: ll.lO.St, 12.38.Lg, 21.65. Qr 

1 Introduction 

Key issues of modern physics of dense matter are concepts explaining the phe- 
nomena related to the appearance of quantum condensates in dense Fermi sys- 
tems. Two regimes are well-known: the Bose-Einstein condensation (BEC) of 
bound states with an even number of Fermions and the condensation of bosonic 
correlations (e.g., Cooper pairs) in the continuum of unbound states according 
to the Bardeen-Cooper-Schrieffer (BCS) theory. While the former mechanism 
concerns states which are well-localized in coordinate space as they occur for 
strong enough attractive coupling, the latter mechanism applies to states which 
are correlated within a shell of the order of the energy gap A around the Fermi 
sphere in momentum space but delocalized in coordinate space. The transi- 
tion between both regimes is called BEC-BCS crossover. Recently, this transi- 
tion regime became accessible to laboratory experiments with ultracold gases 
of fermionic atoms coupled via Feshbach resonances with a strength tunable by 
applying external magnetic fields, see Fig. 1. After the preparation of fermionic 
dimers in 2003, now also the BEC [U [2] and superfluidity of these dimers has 
been observed OH]- The BEC-BCS crossover is physically related 5J to the 
bound state dissociation or Mott-Anderson delocalization transition [6J where 
the modification of the effective coupling strength is caused by electronic screen- 
ing and/or Pauli blocking effects [7]. It is thus a very general effect expected to 
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Figure 1 : Illustration of the transition from BCS pairing with delocalized wave 
functions to BEC of bound states, well-localized in coordinate space, from Ref. 



occur in a wide variety of dense Fermi systems with attractive interactions [H] 
such as electron- hole systems in solid state physics [9 , electron-proton systems 
in the interior of giant planets |10| , deuterons in nuclear matter [11) or diquarks 
in quark matter ^J. The BEC-BCS crossover transition in quark matter is of 
particular theoretical interest due to the additional relativistic regime |13j . 

A systematic treatment of these effects is possible within the path integral 
formulation for finite-temperature quantum field theories. This approach is 
rather general as it is relativistic and is especially suited to take into account the 
effects of spontaneous symmetry breaking. Within these lectures we will present 
the basics of this approach on the example of a model field theory of the Nambu- 
Jona-Lasinio type for a relativistic strongly interacting Fermi system. These 
investigations are also motivated by the analogies of the strongly coupled quark- 
gluon plasma (sQGP) at Relativistic Heavy Ion Collider (RHIC) in Brookhaven 
[15 j with the experiments on BEC of atoms in traps. Furthermore, qualiative 
insights into possible effects observable in the upcoming CBM experiment at 
FAIR Darmstadt as well as from neutron stars with quark matter interiors 
could be derived along the lines of this approach. 



2 Path integral formalism 

2.1 Partition function and model Lagrangian 

As a generic model system for the description of hot, dense Fermi-systems with 
strong, short-range interactions we consider quark matter described by a model 
Lagrangian with four-fermion coupling. The key quantity for the derivation of 
thermodynamic properties is the partition function Z from which all thermo- 
dynamic quantities can be derived. It is given as a path integral which in the 
imaginary time formalism {t — —it) can be expressed as |16j 

Z = jv{iq^)V{q) ef ''"^ i'> , (1) 
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where the chemical potential is introduced as a Lagrange multiplier for assur- 
ing conservation of baryon number as a conserved charge carried by the quarks. 
The notation d^x is shorthand for dr J cPx where (3 = 1/T is the inverse 
temperature. The quark matter is described by a Dirac Lagrangian with internal 
degrees of freedom {Nf = 2 flavors , Nc = 3 colors) , with a current-current-type 
four-fermion interaction inspired by one-gluon exchange 

C = q(i$-mo)q''-—^q—-i^,qq—-l''q, (2) 

a=l 

where A° are the Gell-Mann matrices for color SU{2)). After Fierz transfor- 
mation of the interaction, we select scalar diquark channel and the scalar and 
pseudoscalar channel so that our model Lagrangian assumes the form 

C — Cq + Cqq + Cqq (3) 

where the different terms are given by 

£o q{i$ - mo)q (4) 

C-qq = Gs {qqf + {qilsrqf (5) 
Cqq = GD{q[ij5CT2X2]q^} {q'^ [iC-f5T2X2]q} , (6) 



where 7^ are the Dirac matrices, are SU{2) flavor matrices and C = ij'^^'^ 
is the charge conjugation matrix. Gs and Gd are the coupling strengths corre- 
sponding to the different channels, see Ref. [T7] for a recent review. For the nu- 
merical analysis we adopt parameters from Ref. [18] and consider rjjj = Gd / Gs 
as a free parameter of the model. 

A general method to deal with four-fermion interactions in the path integral 
approach starts with the Hubbard-Stratonovich transformation T9j of the par- 
tition function to its equivalent form in terms of collective bosonic fields, which 
is more suitable to deal with nonperturbative effects such as the occurrence of 
order parameters related to phase transitions in the system as well as collective 
excitations (plasmons = mesons and pairs = diquarks) in these phases. 

2.2 Hubbard-Stratonovich transformation: Bosonization 

The Hubbard-Stratonovich transformation is a two-step procedure which con- 
sists of (1) linearization of the four-fermion interaction terms by introducing 
bosonic auxiliary fields in the appropriate channels and (2) integrating out the 
fermions analytically. 

We introduce the Hubbard-Stratonovich auxifiary fields A{t,x), A*{t,x), 
7r(T, x) and a{T, x) so that the partition function of the system becomes 



Z = / VA*VAVaVTT 



e 



(7) 



dq] [dq] '^'^^{9(^^+f^''fo-mo)q-q(<7+i-f5T-rr)q-^q'^ Rq-^qRq'^) 
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where R — iCj^ T2 (8" A2, R — ijsC (^72(^X2- By introducing Nambu-Gorkov 
spinors 

1 ( q \ = 1 



with q'^{x) = C(f"{x), the Lagrangian takes the bihnear form 

iA*75T2A2 i0 - fijo - m - i^^r ■ tt J 

with 771 — Too + Hence the partition function becomes a Gaussian path 
integral in the bispinor fields which can be evaluated and yields the fermion 
determinant 

Z = j VA*VAVaVTre^ '^^''^^^^ J V^V^el" (10) 

VA*VAVaVTTe^ J "^^^^^ • Det[S'-i] , (11) 

where the inverse bispinor propagator is a matrix in Nambu-Gorkov-, Dirac-, 
color- and flavor space which after Fourier transformation reads 

^-1 ^ f (icj„ + ,u)7o - TO - i7P - «75T • TT ?A75T2A2 \ 

\ iA*j5T2\2 (ZCJ„ - /i)7o - TO - i7p + i75T • TT / 

So far we could derive with (llip a very compact, bosonized form of the quark 
matter partition function ([T]) which is an exact transformation of ([T]), now for- 
mulated in terms of collective, bosonic fields. As we will demonstrate in the 
following, this form is suitable since it allows to obtain nonperturbative results 
already in the lowest orders with respect to an expansion around the stationary 
values of these fields. In performing this expansion, we may factorize the parti- 
tion function into mean field (MF), Gaussian fluctuation (Gauss) and residual 
(res) contributions 

Z{fi,T) = e-'3f^(^'^) - (Ai,T)ZGauss(M,T)Z„s(Ai,T) . 

In the following wc will discuss the physical content of these approximations. 



2.3 Mean-field approximation: order parameters 

In thermodynamical equilibrium, the mean field values satisfy the stationarity 

1 

13V 

dn]\jF dflMF dilniF 



condition of the minimal thermodynamical potential flMF — ^Tfu i.e. 



daMF diTMF dA 



MF 



= , (13) 



equivalent to the fulfillment of the gap equations omf = —^Gs'^t:[Smf) = 
m-mo, TTMF = —HGsTt (jbtSmf) = and Amf = 4G£)Tr {J5T2X2SMF) = A, 
together with the stability criterion that the determinant of the curvature matrix 
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formed by the second derivatives is positive. After the evaluation of the traces 
in the internal spaces and the sum over the Matsubara frequencies one gets 



n 



MF — 



1 [m-moY 



4G 



D 



Tr {In pSjylpf 



(m-mp)^ |A|2 
4Gs 4Gz: 

+ 2rin(l 



(2^ 



S+ + Sp +Sp + 2Tln(l + e- 



p ) + Tln(l + e^^^p ) + rin(l 



(14) 



+ A2 with 

2 The A 7^ dispersion law is associated to the 



where we have defined the particle dispersion relation = \J (^] 

= £'p ± /i, Ep = v/m2 +p 
red and green quarks {E~) and antiquarks (-Ep), whereas the ungapped blue 
quarks (antiquarks) have the dispersion ^p ((Cp ))• From Eqs. (fT3|l with (fT4|) 
we obtain the gap equations for the order parameters m and A, which have to 
be solved self-consistently, 



m — mo 



8Gsm 



(fp 1 



[l-2n;^(i?p)] 



E^ 



[1 - 2np(i?+)] ^ + nf^(-e+) - n;^(ep )| , 



A 



8Gr 



(pp 
(2^ 



2»f(i?p) 



2nF(£;+) 



(15) 
(16) 



with the Fermi dsitribution function np{E) = (1 + e^^) ^ . For zero tempera- 
ture, the gap equations take the simple form 



m — Too 
A 



SGsTO 
8Gd A 



d^p 1 

d^p 
(2^ 



Ei 



Er, 



E^ 



(17) 
(18) 



Solutions of the gap equations for the dynamically generated quark mass to and 
for the diquark pairing gap A at T = as a function of the chemical potential 
are shown in Fig. [51 From the knowledge of the order parameters as functions 
of the thermodynamical variables (T, /x) we can deduce the phase diagram of 
Fig.H 



2.4 Phase diagram 

From the solutions of the gap equations for the order parameters in dependence 
of the thermodynamical variables T and /i we have constructed the phase dia- 
gram of the present quark matter model in the T — fi plane, see Fig. [31 The two 
order parameters allow to distinguish 4 phases: 
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Figure 2: Order parameters for chiral symmetry breaking (full lines) and color 
superconductivity (dashed lines) at T = for different values of the diquark 
coupling rjo- First order phase transitions turn to second order or even crossover 
when riD is increased. For details, see text. 
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Figure 3: Phase diagram of two-flavor quark matter with critical lines for chiral 
symmetry breaking (dashed) and color superconductivity (solid) for three values 
of the diquark couphng strength: r/D = 0.75 (black), 1.0 (green) and 1.3 (red). 
The BEC-BCS crossover occurs when the chiral transition (coincident with the 
Mott transition for mesonic and diquark bound states) occurs inside the 2SC 
phase. It is characterized by the coexistence of diquark condensation with chiral 
symmetry breaking. 
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• A = 0, TO ~ mo: normal phase (NQM) 

• A 7^ 0, TO ~ mo: color superconductor (2SC) 

• A = 0, TO ^ toq: chiral symmetry broken phase (xSB) 

• A 7^ 0, TO 3> Toq: coexistence of xSB and 2SC (BEC phase) 

Order parameters are indicators of phase transitions. The phase transitions can 
be classified according to their order, depending on the behavior of the order 
parameters with the change of thermodynamic variables 

• first order: order parameter jumps, as in the case of the x^B 2SC phase 
transition at not too large coupling. 

• second order: order parameter turns continuously to zero. The 2SC — > 
NQM transition with increasing T is always second order. The xSB 
2SC transition turns from 1st to 2nd order for strong enough coupling. 

• crossover: the order parameter changes continuously, but does not go to 
zero and also has no jumps. An example is the xSB — > NQM transition at 
temperatures above the critical endpoint (CP) where it goes over to the 
line of 1st order transitions in the T — jj, plane. The identification of the 
CP is a key issue for experimental research and thus for the verification 
of QCD models. It is suggested that verifiable signatures (change of the 
fiuctuation spectrum, latent heat or not) are related with it. For strong 
coupling, the CP moves to lower T and finally to T = (for ?7d > 1.3 the 
chiral transition is always crossover). 

Increasing the diquark coupling leads to an increase of the diquark gap 
and therefore a rise in the critical temperature for the second order transition to 
a normal quark matter phase. It shifts also the border between color supercon- 
ductivity (2SC) and chiral symmetry broken phase (xSB) to lower values of the 
chemical potential. For very strong coupling ry/j ~ 1, a coexistence region de- 
velopes, where both order parameters are simultaneously nonvanishing. Under 
these conditions, the phase border is not of first order and therefore no critical 
endpoint can be identified. As we are going to explain in the next section, in the 
xSB phase pion and diquark bound states can exist. At the chiral symmetry 
restoration transition, they merge the continuum of unbound states and turn 
into (resonant) scattering states. When this Mott transition occurs within the 
2SC phase (characterized by a nonvanishing diquark condensate) we speak of 
a BEC-BCS crossover: the condensation of diquark bound states (BEC) turns 
into a condensation of resonances, called Cooper pairs (BCS). 

For lower coupling, the critical point occurs and is shown as a colored dot 
in the phase diagram of Fig. [31 

In the next section we will turn towards the interesting question about the 
quasiparticle excitations in these phases. To this end, we will expand the action 
functional in the partition function up to quadratic (Gaussian) order in the 
mesonic fields and arrive at a tractable approximation for the bosonized quark 
matter model PT|) . 
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2.5 Gaussian fluctuations: bound &: scattering states 

Let us expand now the mesonic fields around their mean field values. In these 
lectures, we will focus on fluctuations in the mesonic channels, where the pion 
and the sigma meson will emerge as quasiparticle degrees of freedom. On the 
example of the pion we will explain the physics of the Mott transition. As 
discussed in the previous section, the phenomenon of the BEC-BCS crossover 
in the 2SC phase is due to the Mott transition for diquarks. The detailed 
investigation of the quantized diquark fluctuations, which are also a prerequisite 
of the formation of baryons, will be given elsewhere [?, [20] . As we already 
noticed, the pion does not contribute to the mean field, and we need to introduce 
only the sigma- field fluctuations as ct — > ctmf + ct- Hence it is possible to 
decompose the inverse propagator into a mean field part and a fluctuation 
part = S^^p + S , where the matrix E is defined as 



(19) 



-(T — 175"? • TT 

—(7 — «75T* • 7? 

In the Gaussian approximation the fermion determinant becomes 
Det \S^^] L ( 1 f d^Q 1 

The propagator Smf is obtained from (fT^ by the matrix inversiorQ 

Smf = ( p+ I ^"^^^ 



with the matrix elements 



tp Pa - tphp po + Sp^p 

n - -EE JLp^ f ,.. Ap-75r.A., (23) 
sp tp ^-^p Po - h^P 

where Sp,ip = ±1, (A+, A~) = (A*, A). For the subsequent evaluation of traces 
in quark-loop diagrams, it is convenient to use this notation with projectors in 
color space, V^g — diag(l, 1, 0),Vh ~ diag(0, 0, 1) and in Dirac space. 



1 ± 70 



7 • p + m 
^P 



The summation over Matsubara frequencies po = icun is most systematic using 
the above decomposition into simple poles in the pq plane. The poles of the 
normal propagators are given by the gapped dispersion relations for the 



paired red-green quarks (antiquarks) , — y (^p )^ -I- A^, and the ungapped 
dispersions — Ep ± /i for the blue quarks (antiquarks). The anomalous 



'^contribution by M. Buballa 
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propagators Fp are only nonvanishing in the 2SC phase when the pair amphtude 
is nonvanishing. Let us notice exphcitly that this procedure has yielded an 
effective action that includes the fluctuation terms responsible for the excitation 
of scalar and pseudoscalar mesonic modes. The evaluation of the traces (^0]) 
can be performed with the result 



-Tr (Smf^Smf^) = (tt, cr) 



n 



TTTT ^ 

n,. 



(24) 



with 



p p+q 



^p p+q 



F+F"j_ 1 

p p+ql 



n,,(go,q) = Tr[G+G++, 
n,,(go,q) = -Tr[G+(75T)G++^(75f)+Fp(75f*)F++,(75f) 
+F+(75f)Fp+^(75f*) + Gp (75f*)Gp+^(75f*)] 



(25) 



(26) 



These polarization functions are the key quantities for the investigation of 
mesonic bound and scattering states in quark matter. In the following we per- 
form the further evaluation and discussion for the pionic modes, the a modes 
is treated in an analogous way. We start with the evaluation of traces and 
Matsubara summation. 



(27r)3 2^ 

Sp,Sk 



npisp^p) - riFiskCp+q) 



90 



p+q 
p+q 



tp,tk %+q 



tptk npitpEp") - npitkEp^^) 



90 



{tptkEp E^j^^ + SpSk^p^]'' 



p+q 



ipEp 



(27) 



where 



T:^{sp,sk) = i 1 



^p^k ~ 



p • (p + q) - ™^ 

EpEp+q 



(28) 



For a pionic mode at rest in the medium (q = 0) this reduces to 



n..(9o,0) = 8 J 



1 - 



E+Ep 



1 + ^ 



E+Ep 



(2^ 

m{e;,e^) 

N{E+E^) 



90 
1 



90 



90 



2E^ 



90 



E- 



90 



Er, 



90 



1 



2Er, 



E^ 



E- 



Er. 



(29) 



where we have introduced the phase space occupation factors N(x, y) = 1 — 
np^x) — np{y) (Pauli blocking) and M{x,y) = np{x) — np{y). For fi this 
function has three poles from the first terms in each bracket, corresponding to 
positive energies (go > 0). So we need to focus only on these three terms. For 
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/U = the second term vanishes due to the prefactor and we are left with two 
poles. 

We make use of the Dirac identity lim^^o ~ ~ ^''"^(x) in order to 
decompose the polarization ftmction into real and imaginary parts after analyt- 
ical continuation to the complex plane. The imaginary part is starightforwardly 
integrated after transformation from momentum to energy oj. At the pole the 
variables transform as 



Pa 



4(^2 _ 4^2^ 



(30) 



For < 1 we know that A = if m > /i, what includes that > 2/i as 
this is the relevant threshold. Therefore, the pole is not hidden and we recover 
the usual 2m threshold. For small enough couplings, A 7^ only if m < /it. 
Therefore, this pole is not hidden in this case. This reasoning includes that the 
argument of the square root is strictly positive. The integration borders thus 
shift p e (0, 00) — > a; e {X±, 00), where the thresholds are given by 2m and 



X± = + m)^ + A2 ± y/{m - /z)2 + A2 . (31) 

The pion polarization function in the 2SC phase can thus be decomposed into 
real and imaginary parts in the following form 



+ 



+ 



1- 



(2^ 



- iV, 0) 
V 



ilmn^^(w + ir?,0) 
1 



2K 



2E^ 



A2 



E+E- 



1 + 



^P^^P- 



E+E^ 



M{E+,E-) 

n{e;,e-) 



V 



E^^E- 
V 



+ E+-E- 



Ep io + E+ + Ep 



-i-{pZEpoN{^+„ , )e(u; - 2m) 



+PwEp^ 



+ -^M{El,E-je{^-X_) 



(32) 



where V denotes the principal value integration, p° = pi^ |a=o= ~ ^^'^ 
we have made explicit the three thresholds, 2m and X±, for the occurrence of 
the corresponding decay processes giving rise to the partial widths and r± , 
respectively. In the normal phase this reduces to 



nOja; + zr?,0) = RenO,(^ 



24 



dp 



(2.)3^(^P'^P) 



V 



2Ep w + 2Ep 



-i-plEpoN{^+ , )e{uj - 2m) 



(33) 
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The analytic properties of the mesonic modes can be analyzed from their spectral 
function. Here we discuss results for pionic modes with q = in the rest frame 
of the medium 



p^(w + i7/, 0) 



(34) 



[1 - 2G,Ren,,(Lj, 0)]2 + [2G^ImU^^{Lj + i7], 0)]2 ' 



In the limit of vanishing imaginary part, we recover the spectral function for a 
"true" (on-shell) bound state 



which corresponds to an infinite lifetime of the state and a mass to be found from 
the pole condition 1 — 2Go-Ren^7r("i7r, 0) — 0. In Fig. [H we show results for 
the mass spectrum of pions and sigma-mesons as a function of the temperature 
for vanishing chemical potential /zs = and strong diquark coupling rjD ~ 
1.0. Since A = 0, the only threshold for the imaginary parts of meson decays 
is 2m. The a mass is always above the threshold and therefore this state is 
unstable in the present model. The pion, however, is a bound state until the 
critical temperature for the Mott transition TMott = 212.7 MeV is reached. For 
T > TMott the pion becomes unstable for decay into quark-antiquark pairs. 
As can be seen from the behavior of the spectral function in the lower panel of 
Fig. in the pion is still a well-identifyable, long-lived resonance in that case. The 
detailed analytic behavior of the pion at the Mott transition has been discussed 
in the context of the NJL model by Hiifner et al. [21], see also the inset of the 
lower panel of Fig. [D It shows strong similarities with the behavior of bound 
states of fermionic atoms in traps when their coupling is tuned by exploiting 
Feshbach resonances in an external magnetic field, see [22] ■ In the context of 
RHIC experiments, one has discussed such quasi-bound states as an explanation 
for the perfect liquid behavior of the sQGP |23j . 

Next we want to discuss the pionic excitations in the presence of a diquark 
condensate in the 2SC phase, see Fig. |3| We choose /i = 320 MeV and discuss 
the effect of melting the 2SC diquark condensate by increasing the temperature 
from T = to T > Tc, where Tc = 95 MeV is the critical temperature for 
the second order transition to the normal quark matter phase. We observe the 
remarkable fact that the 2SC condensate stabilizes the pion at T = as a true 
bound state, although the pion mass exceeds by far the threshold 2m. This effect 
is due to a compensation of gapped and ungapped quark modes and has been 
discussed before by Ebert et al. 20J for T = only. Here we extend this study 
to the finite temperature case, where the pion is still a very good resonance, but 
obtains a finite width. At the critical temperature Tc, the normal pion width is 
restored. But already before T = Tc is reached, the threshold X- is reached and 
the corresponding decay process is opened with a considerable width of 0(50 
MeV). From the pion spectral function in the lower panel of Fig. |3| we observe 
the gap in the excitation spectrum due to the presence of the diquark gap. At 
T > Tc, a resonance type spectral function with a threshold at a; = 2m and a 
resonance peak at w ~ 250 MeV is obtained. 

The discussion of the mesonic modes in the 2SC phase points to a very rich 
spectrum of excitations which eventually lead to specific new observable signals 



lim 



p^{uj + iTj, 0) = 2tt5{1 - 2G<,ReH^^(w, 0)) , 



(35) 
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Figure 4: Upper panel: Mass spectrum of mesons (tt, a) as a function of the 
temperature for vanishing chemical potential fiB — and strong diquark cou- 
pling riD = 1-0. The threshold Eth — 2 uiq for Mott dissociation of pions 
and occurrence of a nonvanishing decay width = Im HT^/m^ is reached 
at Tiviott = 212.7 MeV (see inset). Lower panel: Spectral function for pionic 
correlations for fj,B = in the vacuum at T = (see inset) and at different 
temperatures around the Mott transition. Below TMott, the bound state (delta 
function) and the continuum of scattering states are separated by a mass gap. 
Above TMott , the spectral function is still sharply peaked, related to a lifetime of 
pionic correlations of the order of the lifetime of a fireball in heavy-ion collisions 
(quasi-bound states in the quark plasma). 
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Figure 5: Upper panel: Mass spectrum of mesons (tt, a) as a function of the 
temperature for finite chemical potential /is = 320 MeV and strong diquark 
coupling rjD = 1.0 in the 2SC phase. Below the threshold i?th = 2 niq for 
the onset of the decay width there is another process due to the lower 
threshold i?+ — E- switching on (see inset). Lower panel: Spectral function for 
pionic correlations for /is — 320 in the 2SC phase at T = and at different 
temperatures around the critical temeprature for the 2SC phase. 
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of this hypothetical phase. The CBM experimenid planned at FAIR Darmstadt 
and the NICA projeclH at JINR Dubna could be capable of creating thermody- 
namical conditions for the observation of these excitations in the experiment. In 
view of this discovery potential, we want to outline a few points for the further 
development of the theoretical approach. 

3 Further Developments 

In this contribution we have described the first steps into the interesting and 
very complex physics of the relativistic BEC-BCS crossover theory. As the next 
steps following this development, some of the approximations can be removed. 
In particular, one should next 

• evaluate the full spectrum of diquark states, including their mixing with 
mesonic channels 

• study the backreaction of the correlations on the mcanficld (selfconsistent 
meanfield) 

• include higher orders in the one-fermion-loop approximation (diquark- 
diquark and diquark-meson interactions) 

• study the effect of the color neutrality condition by adjusting color chem- 
ical potential(s) 

• study the effect of charge neutrality (gapless superconductivity) 

• evaluate the contribution from diquark-antidiquark annihilation to the 
photon propagator (Maki-Thompson and Aslamasov-Larkin terms) . 

In particular the latter point bears a big potential for applications to the diag- 
nostics of dense qurk matter formed, e.g., in not too high- energy nucleus- nucleus 
collisions. The onset of color superconductivity not only changes the spectrum 
of diquark states (occurrence of the Goldstone bosons) but due to the nonvan- 
ishing diquark gap additional terms for the diquark annihilation process into the 
observable dilepton channel arise which stem from the nonvanishing anomalous 
propagator contributions. Since the critical temperature for the color supercon- 
ductivity transition might be as high as 100 MeV there is a fair chance to observe 
traces of this transition with the future CBM experiment at FAIR Darmstadt. 

4 Conclusions 

We have derived the gap equations and Bethe-Salpeter equations for the simul- 
taneous treatment of quantized density fluctuations (particle-antiparticle modes 
— mesons) in RPA approximation (one-fermion-loop polarization function) and 
quantized pairing fluctuations (particle-particle modes = diquarks) in ladder ap- 
proximation (no crossed-ladder or vertex corrections) within the path integral 

^contribution by P. Senger 
^contribution by A.S. Serin et al. 
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formalism, i.e. a fully relativistic field theoretical treatment. For the a priori 
unknown interaction, a local current-current coupling (NJL model) has been 
used which in the nonrelativistic limit becomes equivalent to the BCS model 
of superconductivity. After fixing the parameters of the model to the light me- 
son spectrum in the vacuum, the diquark coupling remains as a free parameter 
which has been used to extend the model beyond the traditional range of ap- 
plications into the region of BEC-BCS crossover, where both diquark bound 
states and scattering states occur simultaneously and determine the physical 
properties of the system. Recently, the tuning of the coupling in (nonrelativis- 
tic) low-temperature systems of fermionic atoms with Feshbach resonances in 
an external magnetic field could be used to investigate the BEC-BCS crossover 
in the laboratory. The nonrelativistic limit of the present approach can be used 
to interprete such results within a local coupling model. The fully relativistic 
form can be applied to model strong coupling QCD and to investigate the effects 
of finite temperature and density on the phase structure. We have presented in 
this contribution the phase diagram of quark matter at strong and very strong 
coupling, delineating a possible region of BEC-BCS crossover. The origin of 
the BEC-BCS crossover in superconducting quark matter is the Mott transition 
for diquark bound states. We explain the physics of the Mott transition on 
the example of mesonic correlations. We investigate the spectral function for 
pionic correlations (bound and scattering states) outside (T > Tc) and for the 
first time also inside (T < Tc) the color superconductivity region. We find the 
thresholds for the dissociation of pionic bound states into unbound, but resonant 
scattering states in the quark-antiquark continuum. For the lifetime of pionic 
resonances in the quark matter phase diagram see [H]. Summarizing, we have 
provided the framework for a study of hot and dense fermion matter beyond the 
mean field within a fully relativistic approach, with yet local interaction and in 
Gaussian approximation. The applications to QCD matter within a NJL model 
provide us with a phase diagram for quark matter which is now augmented with 
information about the presence of strong mesonic and diquark correlations with 
possible consequences in the phenomenology of relativistic heavy-ion collisions 
and the interiors of compact stars. 
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